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Magnetic Reconnection with Asymmetry in the Outflow 
Direction 

N. A. Murphy, 1,2,3 C. R. Sovinec, 1,4 and P. A. Cassak 5 

Abstract. Magnetic reconnection with asymmetry in the outflow direction occurs in 
the Earth's magnetotail, coronal mass ejections, flux cancellation events, astrophysical 
disks, spheromak merging experiments, and elsewhere in nature and the laboratory. A 
control volume analysis is performed for the case of steady antiparallel magnetic recon- 
nection with asymmetric downstream pressure, which is used to derive scaling relations 
for the outflow velocity from each side of the current sheet and the reconnection rate. 
Simple relationships for outflow velocity are presented for the incompressible case and 
the case of symmetric downstream pressure but asymmetric downstream density. Asym- 
metry alone is not found to greatly affect the reconnection rate. The flow stagnation point 
and magnetic field null do not coincide in a steady state unless the pressure gradient is 
negligible at the flow stagnation point. 



1. Introduction 

While most two-dimensional models of magnetic recon- 
nection assume that the process is symmetric to a 180° 
rotation about the X-point, there are many situations in 
nature and in the laboratory where this assumption is in- 
valid. In recent years, many papers have addressed mag- 
netic reconnection with asymmetry in the inflow direction 
[e.g., La Belle-Hamer et al, 1995; Nakamura and Scholer, 
2000; Swisdak et al., 2003; 0ieroset et ai, 2004; Borovsky 
and Hesse, 2007; Cassak and Shay, 2007, 2008, 2009; Birn 
et ai, 2008; Murphy and Sovinec, 2008; Mozer et ai, 2008; 
Pritchett, 2008; Borovsky et ai, 2008; Tanaka et al, 2008; 
Mozer and Pritchett, 2009; Eriksson et al., 2009]. In partic- 
ular, Cassak and Shay [2007] generalized the Sweet-Parker 
model [Parker, 1957; Sweet, 1958] to account for reconnec- 
tion between plasmas with different upstream densities and 
magnetic field strengths. They found that the reconnection 
rate is governed by a hybrid Alfven speed that takes into 
account the densities and magnetic field strengths for the 
two upstream regions. The positions of the magnetic field 
null and flow stagnation point are displaced from each other, 
with the field null position set by balance of energy flux and 
the stagnation point position set by balance of mass flux. In 
addition to reconnection with asymmetric inflow, there are 
many situations in nature and the laboratory for which the 
system is asymmetric in the outflow direction. In this paper, 
we analyze steady magnetic reconnection with asymmetry in 
the outflow direction. 

The best known scenario for magnetic reconnection with 
asymmetry in the outflow direction is the Earth's magne- 
totail. In this case, asymmetry is a particularly important 
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consideration because it helps determine the amount of en- 
ergy transported in the earthward and tailward directions 
as a result of reconnection. At distances of ~5-15_Re, there 
is a considerable pressure gradient as the plasma pressure 
decreases approximately monotonically with distance from 
Earth [Lui et ai, 1994; Shiokawa et al., 1997; Xing et al., 
2009]. Earthward-directed reconnection outflow must work 
against strong gradients in both plasma pressure and mag- 
netic pressure. Because of the global configuration of the 
magnetotail, the X-line characteristically moves in the tail- 
ward direction [Hones, 1979]. Reconnection with asymme- 
try in the outflow direction has often been seen in simula- 
tions of the magnetotail [e.g., Birn et al., 1996; Hesse et al., 
1996; Hesse and Schindler, 2001; Kuznetsova et al., 2007; 
Laitinen et ai, 2005; Laitinen, 2007; Birn and Hesse, 2009; 
Zhu et al., 2009], though the degree of asymmetry depends 
on the proximity of the reconnection layer to Earth and how 
reconnection is driven. The largest discrepancy between 
earthward and tailward outflow velocities in these simula- 
tions was seen by Laitinen et al. [2005] and Laitinen [2007], 
where the inflow had a large component of velocity in the 
outflow direction; consequently, there is a large separation 
between the X-line and the flow reversal line in their results. 

Observing reconnection with asymmetry in the outflow 
direction in the magnetosphere requires multiple satellites 
crossing the earthward and tailward sides of the diffusion 
region at approximately the same time. While statistical 
approaches are possible [e.g., Petrukovich et al., 2009], ob- 
servations of a single event are not common. One occurrence 
is a crossing of the diffusion region by Cluster on 11 Octo- 
ber 2001. Cluster was in the region between the outflow 
jets from 03:30-03:36 UT, but passed the X-line at 03:31 
UT. One possible explanation of this is that the X-line was 
near the tailward end of the diffusion region. However, other 
explanations (e.g., time-dependent behavior or undetected 
additional X-lines) cannot be ruled out with the available 
data [Laitinen et ai, 2007]. 

In solar physics, reconnection during coronal mass ejec- 
tions (CMEs), solar flares, and flux cancellation events are 
asymmetric in the outflow direction when one outflow jet 
is directed sunward and the other outflow jet is directed 
away from the Sun [e.g., Kopp and Pneuman, 1976; Martin 
et al, 1985; Shibata et al., 1995; Litvinenko, 1999; Lin and 
Forbes, 2000; Aurass et ai, 2002]. Observations of bidirec- 
tional jets in the solar atmosphere [see Innes et al., 1997; 
Wang et al, 2007; Liu et al, 2009; Gontikakis et al, 2009] 
show that, despite the effects of gravity, the redshifted jet 
is often slower than the blueshifted jet because the red- 
shifted jet must propagate into a higher density medium. 
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In these events, gravity's most important effect is the es- 
tablishment of a stratified medium. Current sheets forming 
in such a medium are likely to have strong gradients in the 
outflow direction for upstream density, pressure, and mag- 
netic field strength [see Ciaravella et al, 2002; Ko et al, 
2003; Chen et al, 2004; Bemporad et al, 2006; Lin et al, 
2007; Ciaravella and Raymond, 2008; Bemporad, 2008; Lin 
et al., 2009; Vrsnak et al, 2009; Saint-Hilaire et al., 2009; 
Aurass et al, 2009]. Simulations of reconnection in a strat- 
ified medium show that the redshifted jet can be up to an 
order of magnitude slower than the blueshifted jet [Roussev 
et al., 2001], and that reconnection in such an atmosphere 
displays a more complicated velocity structure than sym- 
metric two-dimensional reconnection [Galsgaard and Rous- 
sev, 2002]. Gravity itself can be an important consideration 
if the work done by electromagnetic forces is comparable to 
or less than the work done against gravity [Reeves, 2006]. 
Asymmetry in the outflow direction also happens when mag- 
netic field lines in one downstream region are line-tied while 
magnetic field lines in the other downstream region are open. 

During turbulent reconnection [e.g., Lazarian and Vish- 
niac, 1999] and reconnection occurring during a turbulent 
cascade [e.g., Servidio et al., 2009], there will in general be 
many reconnection sites throughout the volume of interest. 
Reconnection occurring at each of these sites will in gen- 
eral be asymmetric in the inflow and outflow directions, as 
well as the out-of-plane direction. Reconnection processes 
involving multiple competing reconnection sites or multiple 
magnetic islands [e.g., Lee and Fu, 1986; Drake et al, 2006; 
Lin et al, 2008a; Chen et al, 2009] will also likely involve 
asymmetry in the outflow direction, especially if the X-lines 
are not evenly spaced or develop at different rates. 

In astrophysical settings, the winds of strongly magne- 
tized hot stars (e.g., the Bp star a Ori E) can be channeled 
along a predominantly dipolar field to form an equatorial 
circumstcllar disk or buildup of material [Nakajima, 1985; 
Cassinelli et al, 2002; Townsend and Owocki, 2005]. While 
the dipole field is in general dominant close to the star, re- 
cent axisymmetric simulations show that the continual tun- 
neling of material can eventually lead to centrifugal break- 
out events associated with magnetic reconnection [ud-Doula 
et al, 2006, 2008]. In this case, the reconnection outflow 
is aligned with the radial direction, with one exhaust path 
directed towards the disk and the star, and the other lead- 
ing to the interstellar medium. Such reconnection events 
could be the source of the X-ray flares observed on a Ori E 
by ROSAT [Groote and Schmitt, 2004]. Considerations of 
asymmetry in the outflow direction are also important for 
magnetic reconnection events associated with centrifugal in- 
stabilities and plasma release in the Jovian magnetosphere 
[e.g., Kivelson and South-wood, 2005]. 

In the laboratory, reconnection with asymmetry in the 
outflow direction occurs during the merging of sphcromaks 
and in toroidal plasma configurations where the reconnec- 
tion outflow is aligned with the radial direction. Relevant ex- 
periments include the Swarthmore Spheromak Experiment 
(SSX) [Cothran et al, 2003], the Magnetic Reconnection 
Experiment (MRX) [Yamada et al, 1997], and TS-3/4 at 
the University of Tokyo [Ono et al, 1993]. Recent sphero- 
mak merging experiments at MRX have shown that asym- 
metry in the outflow direction develops as a result of the 
Hall effect [Inomoto et al, 2006]. In these experiments at 
MRX, the reconnecting magnetic field lines do not lie in 
the poloidal plane, and there is a component of the elec- 
tron flow associated with the reconnection current in the 
radial direction. This radial component of electron velocity 
pulls the reconnecting field lines, leading to a shift in po- 
sition of the X-point, asymmetric outflow, and asymmetric 
downstream pressure. Reversing the toroidal field direction 
changes the direction of the shift, but because of toroidicity, 
this also changes the reconnection rate and radial pressure 



profile [Inomoto et al, 2006; Murphy and Sovinec, 2008]. 
Recent simulations of spheromak merging in SSX show re- 
connection with much stronger radially inward directed out- 
flow even though the plasma pressure near R — is large 
due to a pileup of exhaust [Lin et al, 2008b]. These results 
suggest that considerations of asymmetry in the outflow di- 
rection are important for the interpretation of bidirectional 
jets recently reported in experiment [Brown et al, 2006]. 

Murphy and Sovinec [2008] presented simulations of the 
reconnection process in the geometry of MRX, showing that 
asymmetric inflow occurs during the pull mode of operation 
and asymmetric outflow during the push mode of operation 
[see Yamada et al, 1997, Figure 3]. The inboard (low radius) 
side of the current sheet is more susceptible to buildup or 
depletion of density due to the lesser available volume than 
on the outboard (high radius) side of the current sheet. As a 
result of the pressure buildup at low radii during push recon- 
nection, the X-point is located closer to the outboard side of 
the current sheet than the inboard side. Consequently, the 
radially inward directed outflow is subjected to a stronger 
tension force than the radially outward directed outflow, al- 
lowing comparable outflow velocities from both the inboard 
and outboard sides of the current sheet (a similar effect is 
discussed by Galsgaard and Roussev [2002]). During several 
time intervals in these simulations and despite the higher 
pressure in the inboard downstream region, the radially in- 
ward directed outflow speed is found to be greater than the 
radially outward directed outflow speed. The magnetic field 
null and flow stagnation point are separated during both 
pull and push reconnection [Murphy, 2009, Figures 2.4 and 
2.6]. Push reconnection is an example of how asymmetry in 
the outflow direction develops when outflow in one down- 
stream region is confined more effectively than outflow in 
the other downstream region. 

Oka et al. [2008] performed particle-in-cell (PIC) simu- 
lations of reconnection where outflow from one end of the 
current sheet is impeded by a hard wall while outflow from 
the other end encounters no such obstruction. They found 
that the X-line retreats from the wall at ~10% of the up- 
stream Alfven velocity Va and that the reconnection rate 
is largely unchanged from the symmetric case. Moreover, 
there is a separation between the ion flow stagnation point 
and the magnetic field null, with the field null located fur- 
ther from the wall than the ion flow stagnation point. In 
a work that relates asymmetry in the inflow direction with 
asymmetry in the outflow direction, Swisdak et al. [2003] 
found that the presence of a density gradient in the inflow 
direction across a current sheet can lead to a drift of the 
X-line in the electron diamagnetic drift direction when a 
guide field is present [see also Rogers and Zakharov, 1995]. 
The reconnection process is suppressed when the drift ve- 
locity is comparable to or greater than the Alfven velocity. 
The effects of current sheet motion and time-dependence 
on slow shock mediated reconnection layers have also been 
considered [Owen and Cowley, 1987a, b; Kiehas et al, 2007; 
Kiehas et al, 2009]. 

In this paper, we perform a control volume analysis for 
a current sheet with asymmetric downstream pressure and 
test the resulting scaling relations against simulations. The 
objectives are to determine (1) the relationship between the 
upstream parameters, the downstream pressures, and the re- 
connection outflow velocity, (2) how the reconnection rate is 
affected by asymmetric downstream pressure, and (3) what 
sets the positions of the magnetic field null and flow stagna- 
tion point. In section 2, we write the equations of resistive 
magnetohydrodynamics (MHD) in a time-independent inte- 
gral form that is amenable to a control volume analysis. In 
section 3, we review the effects of symmetric downstream 
pressure on antiparallel reconnection and develop scaling 
relations for a current sheet with asymmetric downstream 
pressure. In section 4, we test the scaling relations derived in 
section 3 against resistive MHD simulations of reconnection 
with asymmetry in the outflow direction. In section 5, wc 
provide a discussion and summarize our results. A similar 
analysis for reconnection in cylindrical geometry with out- 
flow aligned with the radial direction is presented by Murphy 
[2009, section 3.4]. 
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2. Equations of Magnetohydrodynamics 

The equations of resistive MHD in conservative form [e.g., 
Goedbloed and Poedts, 2004, pp. 165-166] are 
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where B is the magnetic field, E is the electric field, V is 
the bulk plasma velocity, J is the current density, p is the 
plasma pressure, p is mass density, rj is the plasma resistiv- 
ity, £ = pV' 2 /2 + p/(-y - 1) + B 2 /2p is the total energy 
density, and 7 is the ratio of specific heats. The identity 
dyadic tensor is given by I = xx + yy + zz. Equation (3) 
includes the the internal energy flux, pV/(7 — 1), and the 
mechanical work done on or by the plasma by pressure gra- 
dients while moving, pV. 

Following the approach presented by Cassak and Shay 
[2007], we assume a steady-state system, integrate over an 
arbitrary closed volume V bounded by the surface 5, and use 
the divergence theorem to write the continuity, momentum, 
and energy equations as 
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a current sheet with asymmetric downstream pressure while 
relaxing their assumptions regarding compressibility [see 
also Parker, 1963; Chae et ai, 2003; Litvinenko and Chae, 
2009]. The characteristic parameters used in this derivation 
are: Bi n , upstream magnetic field strength; Vi„, plasma in- 
flow velocity; V ut, plasma outflow velocity; p; n , upstream 
plasma pressure; p ut, downstream plasma pressure; pj„, up- 
stream plasma density; pout, downstream plasma density; 
J y , out-of-plane current density inside the layer; E y , out- 
of-plane electric field; L, current sheet half-length; and S, 
current sheet half-thickness. We define x as the outflow di- 
rection, y as the out-of-plane direction, and z as the inflow 
direction. 

Everywhere except within the reconnection layer, the 
ideal Ohm's law is approximately valid. By assuming a 
steady state the electric field is constant and given by 



Ey — VinB. 



(11) 



Since B x reverses over a distance of ~2<5, Ampere's law gives 



V ~ Mo^' 



(12) 



Matching the resistive electric field inside the layer with the 
ideal electric field outside the layer gives 



JL 

p S' 



(13) 



Evaluating the conservation of mass relation given in 
equation (7) over the entire volume of the current sheet 
yields the relation 



pinVi„L ~ poutVoutS. 



(14) 



The conservation of momentum surface integral given in 
equation (8) is satisfied by any distribution of fluxes with the 
assumed symmetry when integrating over the outer bound- 
ary of the current sheet. Evaluating the conservation of 
energy relation given in equation (9) yields the relation 



where dS is a differential area element pointing in the out- 
ward normal direction to S. Similarly, with the help of 
Stokes' theorem, equation (4) leads to 



dS x E = 0. 



(10) 



Equations (7)-(10) are valid for any arbitrary closed volume, 
provided a steady-state has been achieved. These surface 
integrals are evaluated in section 3 to investigate magnetic 
reconnection with asymmetry in the outflow direction. 
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where a = 7/(7 — 1). Here we neglect contributions from 
upstream kinetic energy and downstream magnetic energy. 
Dividing equation (15) by equation (14) and rearranging 
gives the scaling relation 



V 2 



vl 



Pout 
pout 



Pn 
Pi; 



(16) 



3. Scaling Relations 

The Sweet-Parker model [Sweet, 1958; Parker, 1957] de- 
scribes symmetric steady-state antiparallel magnetic recon- 
nection in the resistive MHD framework when compressibil- 
ity, viscosity, and downstream pressure are unimportant. In 
this section, we extend these results to account for reconnec- 
tion with asymmetric downstream pressure. After reviewing 
the effects of symmetric downstream pressure on the recon- 
nection process in subsection 3.1, we consider the case of 
asymmetric downstream plasma pressure in subsection 3.2. 
We then investigate the internal structure of such a current 
sheet in subsection 3.3. 



Pin Mn 



B in Pir 



Pl Pl 








Pr Pr 








— Lr — 


V r 



Pin Mn 



B in Pir 



2L 



25 
\ 



3.1. Effects of Symmetric Downstream Pressure 

The effects of symmetric downstream pressure on a 
Sweet-Parker current sheet are discussed by Priest and 
Forbes [2000, pp. 123-126]. Presently, we review their re- 
sults using the approach that we employ later this section for 



Figure 1. Sweet-Parker-like reconnection with asym- 
metric downstream pressure and a pressure gradient 
through the current sheet. The solid vertical bar inside 
the current sheet (marked s) represents the flow stagna- 
tion point and the dashed vertical bar (marked n) repre- 
sents the magnetic field null. 
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where Va = -Bin /^/po pin is the upstream Alfven speed and 
we ignore factors of order unity. The term ap/ p is the en- 
thalpy per unit mass. Using equation (13), the scaling for 
the dimensionless reconnection rate can then be written as 
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where S = poLVa/v is the Lundquist number. The recon- 
nection rate depends weakly on the downstream pressure 
except when the bracketed quantity is close to zero. The 
Sweet-Parker scalings of V ou t ~ Va and Vi„/VA ~ 5 -1 ^ 2 
are recovered when pout/pin and the quantity in brackets 
are independent of S. We also see that when compression 
makes the outflow density larger than the inflow density, it 
relaxes the usual bottleneck from flow moving through the 
reconnection region. 

3.2. Effects of Asymmetric Downstream Pressure 

We now consider a current sheet with symmetric inflow 
but with asymmetric outflow and downstream pressure. In 
this framework, it is necessary to assume that the current 
sheet position and structure is steady within the inertial ref- 
erence frame of the X-line. For example, reconnection could 
be externally driven in such a way that constrains the po- 
sition of the current sheet. The setup of this problem is 
shown in Figure 1. Throughout this analysis, subscripts L 
and R indicate that the variable represents the characteris- 
tic downstream value of a field for the left and right sides of 
the current sheet. 

To proceed, we evaluate the surface integrals given in 
equations (7)-(9) over the whole volume of the current sheet 
depicted in Figure 1. The conservation of mass surface in- 
tegral given in equation (7) yields the relation 



2Lp in V in ~ plVlS + prVrS. 



(18) 



Evaluating the component of the conservation of momen- 
tum surface integral given in equation (8) in the outflow 
direction yields a relation between the plasma pressures and 
momentum fluxes from each exit of the reconnection layer, 



PlV l + p L ~ p R V R + p R . 



(19) 



Because the current sheet is assumed to be long and thin, the 
above relation neglects forces due to the downstream mag- 
netic field. However, magnetic tension does not need to be 
negligible throughout the volume of integration for this re- 
lationship to hold. Rather, tension need only either be neg- 
ligible along the boundary or contribute along the boundary 
evenly in both outflow directions. If the upstream magnetic 
field is not parallel to the boundaries along z — ±5 in a way 
which is not symmetric in the outflow direction, this may 
yield an additional contribution by tension towards momen- 
tum balance in the outflow direction. Downstream mag- 
netic pressure can be important when the global magnetic 
field configuration contains a large vertical component that 
impedes outflow from one side of the current sheet [e.g., 
Inomoto et ai, 2006, Figure 5]. We also assume that the 
momentum flux pVV into the current sheet does not sig- 
nificantly contribute to momentum balance in the outflow 
direction; this is expected to be important only when the 
outflow component of the inflow velocity is of the same or- 
der as the outflow velocities. Force balance must be met 
in both the inflow and outflow directions simultaneously in 
order for the assumption of time- independence to be valid. 

Using the expression for the electric field given in equa- 
tion (11), the energy conservation integral (9) provides the 
relation 
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The above relation neglects upstream kinetic energy and the 
Poynting flux out of the layer. 

By using equation (18) to eliminate 2LU;„ from equation 
(20) and equation (19) to eliminate Vr, we arrive at the 
following cubic relationship which can be solved for V£ , 



~ C 6L V£ + Ca L v£ + C 2L Vl + Col, 



(21) 



where we do not explicitly assume the nature of the dissipa- 
tion mechanism. The coefficients for the above equation are 
functions of the upstream magnetic field strength as well as 
the upstream and downstream densities and pressures, and 
are given by 
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and we define the average downstream pressure p and the 
pressure difference Ap as 



P = 



PL + PR 



Ap = PR - p L • 



(27) 
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Equation (21) was derived assuming that the scaling fac- 
tors given in equations (18), (19), and (20) are unity. If 
this is not the case, then if £ is equal to the right hand 
side divided by the left hand side of equation (18), and £ is 
equal to the right hand side divided by the left hand side of 
equation (20), then the transformation cf„ — > (C/£) cf„ will 
algebraically account for scaling factors that are not unity 
in equations (18) and (20) for equation (21). 

Equation (21) simplifies for some special cases. When 
Pl = Pr = Pout, the coefficient C%l vanishes, leaving a 
quadratic equation in V^. In the incompressible limit with 
pin — pL — Pr = p and a = 1, the solution becomes 
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Ap 
~2p~ 
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(29) 



where the plus and minus signs refer to Vl and Vr, respec- 
tively. This gives the expected result that the outflow speed 
is slower on the side with higher downstream pressure. 

Next, consider the special case with pl = pr = jw, but 
where the downstream densities can be different. In this 
case, Col vanishes, again leaving a quadratic equation. The 
solution is 
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* , 
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(30) 



where we ignore factors of order unity. This equation shows 
that the outflow speed is higher on the low density side. The 
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Figure 2. Contours of the outflow velocity from the left 
side of the current sheet, Vl/Va, calculated from equa- 
tion (29) as a function of Pl — Pin and p^ — Pin for the 
incompressible case. Contours are separated by 0.2. Fig- 
ures 2, 3, and 4 assume that the scaling factors for equa- 
tions (18), (19), and (20) are unity. 



scalings presented in equations (29) and (30) both reduce to 
the scaling in equation (16) in the symmetric limit as they 
should. 

Equations (21), (29), and (30) are derived solely from 
the scaling relations for conservation of mass, energy, and 
momentum without explicitly assuming a dissipation mech- 
anism [see also Cassak and Shay, 2007]. By using equations 
(13) and (18) (and consequently assuming uniform resistive 
dissipation inside the current sheet), the inflow speed can 
be written as 

vi-y^vis- (31) 

Using equations (11) and (31), the reconnection rate is 
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V (Vl + V R ) 



2fi L 
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Solutions of equation (29) for Vl as a function of pl — pin 
and pr — p in in units of the upstream magnetic pressure are 
presented in Figure 2 for the incompressible limit. The value 
for Vr can be found by switching the values for pl — Pin and 
Pr — Pin- We see that the outflow velocity from one end 
does not depend strongly on the downstream pressure from 
the opposite end. In fact, reconnection events (e.g., in the 
solar atmosphere) do not require bidirectional outflow jets 
traveling at the Alfven speed. Rather, in the presence of 
asymmetric downstream pressure, there can be one Alfvenic 
jet and one sub- Alfvenic jet [see also Roussev et al., 2001]. 
As shown by the widely spaced contours in F igure 3, the nor- 
malized reconnection rate S 1 ^ 2 V in /VA = yj (Vl + Vr)/2Va 
is only weakly dependent on the difference in downstream 
pressures. This conclusion is consistent with the simulations 
of X-line retreat reported by Oka et al. [2008], in which the 
reconnection rate is not greatly affected when outflow from 
one reconnection jet is impeded by the presence of an ob- 
stacle while the other outflow jet has no such obstruction. 



Figure 3. Contours of the normalized reconnection rate, 
given by S 1/2 Vm/V A = V( v l + Vr)/2Va, as a function 
of pl — Pin and pr — pi n and calculated using equation 
(29) to find Vl and Vr for the incompressible case. Con- 
tours are separated by 0.2. 



Figure 4 shows solutions for the incompressible case as 
a function of pl — Pin for fixed PR Pin — Bi/2^o. The 
outflow velocities, calculated using equation (29) and shown 
in Figures 4a and 4b, illustrate the weak dependence that 
downstream pressure from one side of the current sheet has 
on the outflow velocity from the other side of the current 
sheet. Figures 4c and 4d consider the limiting case where 
L is prescribed by external influences on geometry and Vi n 
varies as a function of pl and pr. The normalized reconnec- 
tion rate, seen in Figure 4c, change s modestly despite the 
large change in pl- Defining So = \J r]L/ [ioVa, we see that 
S/S = sJ2Va/{V l + Vr) increases with pl- The increased 
current sheet thickness slows the reconnection rate slightly 
by equation (13). Figure 4e considers a different limiting 
case for which Vi n is prescribed due to external driving of 
reconnection and L varies as a function of pl and pa to 
maintain the same reconnection rate. For this case, S is 
given by equation (13) and is independent of downstream 
pressure. Defining Lq = t)Va/ ^oV 2 n , the normalized length 
is given by L/Lq = (Vl + Vr)/2Va- Figure 4e shows that 
greater downstream pressure reduces the length of the cur- 
rent sheet for this case, and hence the throughput of mass, 
in response to greater downstream pressure. 

3.3. Internal Structure 

Now that the global quantities associated with a current 
sheet with asymmetric downstream pressure can be found, 
we turn our attention to the internal structure of the cur- 
rent sheet. The current sheet is split into three regions of 
lengths Ll, Lm, and Lr, with boundaries at the flow stag- 
nation point and magnetic field null as indicated in Figure 
1. The length Lm is the distance between the magnetic field 
null and the flow stagnation point, which we will see need 
not be zero. The full length of the current sheet, 2L, is given 
by 

2L = Ll + L M + L R . (33) 
We assume in this section that the fields along each bound- 
ary are describable by approximately uniform values for the 
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Figure 4. A solution slice along pr — pin = Bj 2 n /2/io for 
different values of pl — Pin for the incompressible case. 
S hown are (a) V L /V A , (b) V R /V A , (c) S 1/2 V tn /V A = 
V(Vl + Vr)/2Va, (d) 5 /So, and (e) L/L . Plots (c) and 

(d) describe the limiting case where L is prescribed by 
the global geometry and that V% n is free to vary, whereas 

(e) assumes that Vi„ is prescribed externally and that L 
is free to vary. 



upstream fields; however, this may not be justified when 
current sheet motion relative to the upstream fields is im- 
portant or when a long current sheet develops in a stratified 
medium such as the wake behind a CME. 

As in the model by Cassak and Shay [2007] , the position 
of the flow stagnation point is set by conservation of mass. 
Evaluating equation (7) for the three sections of the current 
sheet presented in Figure 1 yields the conservation of mass 
relations 



PinVinLL ~ PlVlS, 
PinVi„LM ~ PnV„S, 
pinVin (Lm + Lr) ~ PrVrS. 



(34) 
(35) 
(36) 



where p n is the density and V n is the outflow component of 
velocity at the magnetic field null. The location of the flow 
stagnation point can be derived from equations (33), (34), 



and (36), and is given by the relations 
PlVl 



Ll 
Lm + Lr 



2L 
' 2L 



PlVl + PrVr 

PrVr 
PlVl + PrVr 



(37) 
(38) 



Evaluating the conservation of energy surface integral 
given in equation (9) in a similar way yields 
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(39) 



(40) 



(41) 



where B 3 is the vertical magnetic field strength at the flow 
stagnation point. When the magnetic field null and flow 
stagnation point are separated, there is a Poynting flux 
across the flow stagnation point and a kinetic energy flux 
across the magnetic field null. 

While, in principle, equations (39)-(41) can be solved for 
Ll, Lm, and Lr, we proceed using an alternate argument to 
find the separation between the magnetic field null and flow 
stagnation point. In a steady state, the outflow component 
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Figure 5. A slice in the outflow direction along z — of 
terms from (a) momentum balance and (b) electric field 
balance for case B at t — 13.3 ps. Vertical bars represent 
the positions of the flow stagnation point (marked 's') 
and the magnetic field null (marked 'n'). 
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of the momentum equation along z 

8V X B z dB x 



pVx- 



dx 



Mo dz 



reduces to 

dp 
dx ' 



(42) 



where we assume that the magnetic pressure inside the cur- 
rent sheet is small compared to magnetic tension. At the 
flow stagnation point, the magnetic tension force must can- 
cel the pressure gradient force in a steady-state system. 
Moreover, the magnetic field null is colocated with the flow 
stagnation point in a steady state only when there is no 
pressure gradient at the magnetic field null. 

To quantify this, define x n as the position of the mag- 
netic field null and x a as the position of the flow stagnation 
point. By definition, x 3 is given by V x (x s ) = 0. Evaluating 
equation (42) at the flow stagnation point using equation 
(12) gives 

B z (x s ) B in ^ dp 
Mo 5 dx 



(43) 



Using that B z (x n 
yields 



0, a Taylor expansion to first order 



B z (x s ) ~ (a 



, dB z 

l) dx 



(44) 



Hence, we can approximate the distance between the flow 
stagnation point and the magnetic field null for a steady 
state, 



(fioS\ ( - d P /dx\ x=1 
\B in )[ 8B z /dx\ x=x 



(45) 



4. Comparison to Simulations 

In this section, we test the scaling relations and results 
derived in Section 3 using resistive MHD simulations for the 
configuration of MRX, but with the geometry straightened 
from cylindrical to linear. Asymmetry in the outflow direc- 
tion develops for the push mode of operation [see Yamada 
et al., 1997, Figure 3] when one downstream wall is closer 
to the current sheet than the other downstream wall. The 
driving process in these simulations constrains the position 
of the reconnection layer and thus limits time-dependent 
motion of the current sheet. Except as otherwise noted, 
the numerical method and simulation setup are described 
by Murphy and Sovinec [2008]. 



4.1. Numerical Method and Problem Setup 

The NIMROD code (Non-Ideal Magnetohydrodynamics 
with Rotation, Open Discussion) [Sovinec et at, 2004] has 
been successfully used to model reconnection in a variety of 
geometries [e.g., Hooper et al., 2005; Murphy and Sovinec, 
2008; Zhu et al, 2009]. NIMROD uses a finite element ex- 
pansion in the poloidal plane and, in three-dimensional sim- 
ulations, a Fourier representation for the out-of-plane direc- 
tion. The system of equations solved by NIMROD for the 
two-dimensional simulations reported in this section are 



^ = -Vx(i,J-VxB)+ KMvtV (V ■ B) , , 
at 

/ioJ = V x B, 

p (j^- + V ■ VV^j = JxB-Vp + V- puW, 
— + V ■ (nV) = V ■ DVn, 

^(§ + v.vr)=-fv.v-v. q + Q, 



(46) 
(47) 
(48) 

(49) 
(50) 



where the heat source term Q — r)J 2 + ^pVV T :VV in- 
cludes Ohmic and viscous heating, isotropic thermal con- 
duction q = — nxVT is used, v is the kinematic viscosity, 
X is the thermal diffusivity, and D is an artificial number 
density diffusivity. Divergence cleaning is included in equa- 
tion (46) to control divergence error [Sovinec et al., 2004]. 
The ratio of specific heats is given by 7 = 5/3, which cor- 
responds to a = 7/(7 — 1) = 2.5. The resistivity is given 
by ri/fio — 20 m 2 s _1 , which corresponds to a Lundquist 
number of S ~ 100 using the half-length of the current 
sheet and the immediately upstream magnetic field strength 
and density. The magnetic Prandtl number is given by 
Pm = v/(ri/fio) = 0.25, with i/ = D = x = 5m 2 s" 1 . 
Reconnection is driven by applying an electric field on the 
surfaces of two flux cores which have a minor radius of 9.4 cm 
and are located at (x, z) = (0 cm, ±27 cm) using the coor- 
dinate system established in section 3. In linear geometry, 
the flux cores are infinite cylinders. The upper and lower 
boundaries are at z = ±62 cm. For one run, the down- 
stream boundaries are at x = —22.5 cm and x = 32.5 cm 
(case A). For the another run, the downstream boundaries 
are at x — —17.5 cm and x = 32.5 cm (case B). 
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Figure 7. Comparisons between model predictions and simulation results. Shown are (a) the outflow 
velocity Vl compared against equation (21), (b) the outflow velocity Vr compared against equation (21), 
and (c) the electric field E y compared against equation (32). The symbol usage is described in Figure 6. 



4.2. Simulation Results 

Next, we present results from these simulations and com- 
pare them to the model presented in this paper. A cut along 
z = from case B at 13.3 /xs is shown in Figure 5. At this 
time, the flow stagnation point is at x = 1.42 cm and the 
magnetic field null is at x = 1.82 cm, indicating a short 
separation between the two points. Over most of the sim- 
ulated time the flow stagnation point is located closer to 
the side with the impeded outflow than the magnetic field 
null in qualitative agreement with equation (45). Magnetic 
pressure is not important within the current sheet. 

To perform quantitative comparisons with theory, the rel- 
evant quantities must be extracted from the numerical re- 
sults. The full length 2L of the current sheet is taken to be 
the distance along z = between the two locations where 
the out-of-plane current density drops to a fraction / of its 
peak value, where / is either 1/e or 1/e 2 . The thickness 
of the current sheet, S, is taken to be the distance in the z 
direction between the location where the out-of-plane cur- 
rent density peaks and where it falls off to / of its peak 
value. The values for the upstream fields are extracted from 
the simulation at z = ±r5 above and below where the cur- 
rent density peaks. This method slightly but systematically 
underestimates the upstream magnetic field strength. The 
values for the downstream fields are taken at z — where 
the out-of-plane current density falls to / of its peak value. 

Comparisons between simulation and our scaling relations 
are shown in Figure 6, using both of the aforementioned val- 
ues of / for both case A and case B. Figure 6a compares the 
left and right hand sides of equation (18) which approx- 
imates conservation of mass, Figure 6b compares the left 
and right hand sides of equation (19) regarding momentum 
balance, and Figure 6c compares the left and right hand 
sides of equation (20) which approximates conservation of 
energy. Verification of these scaling relations requires that 
the data reasonably fit a straight line through the origin. In 
Figures 6a, 6b, and 6c, we see that the left and right hand 
sides of each of the equations approximately fit straight lines 
through the origin with slopes close to unity. In Figures 6a 
and 6c, the slope is slightly greater than unity (~ 1.15-1. 2). 

Comparisons between the outflow velocities extracted 
from simulation and calculated as roots of equation (21) 
are shown in Figures 7a and 7b. Because the positions (and 
even existence within the set of real numbers) of roots of 
high order polynomials can be sensitive to small changes in 
the coefficients [i.e., Wtlkinson, 1959], modest differences be- 
tween the left and right hand sides of equations (18), (19), 
or (20) sometimes lead to large errors in the solution for 



Vl and Vr or the relevant root becoming complex. Be- 
cause of this property common among high order polyno- 
mials, not all of the instances considered have real roots 
and there is increased scatter in Figures 7a and 7b beyond 
what is seen in Figure 6. Despite this, the simulation results 
show reasonable agreement for instances where the roots of 
the polynomials are not greatly impacted by scaling factors 
that are not unity in equations (18), (19), and (20). The 
expression for the reconnection electric field strength given 
by equation (32) is compared against simulation in Figure 
7c, showing good agreement despite a small underpredic- 
tion of ~10-20%. The positions of the flow stagnation point 
given by equations (37) and (38) are tested against simula- 
tion in Figure 8. Despite some outliers, most of the data 
points show a good correspondence between the model pre- 
dictions and the simulation results. The scatter is primarily 
due to time-dependent effects and the non-uniformity of the 
upstream fields. However, the presence of a local pressure 
maximum near the flow stagnation point and magnetic field 
null complicates the determination of Lm for most cases 
because the pressure gradient varies significantly in this re- 
gion; consequently, equation (45) does not reliably predict 
the separation between the flow stagnation point and mag- 
netic field null for these cases. 

As a final check for the assumptions of this model, the 
surface integrals in equations (7), (8), and (9) are calculated 
along the current sheet boundaries using the finite element 
basis functions to interpolate the data. Evaluating the con- 
servation of mass integral in equation (7) shows that the 
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served in simulation compared against the relation given 
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mass influx is within ~ 10-25% of the mass efflux, indicat- 
ing modest time-dependence. Evaluating the conservation 
of energy integral given in equation (9) shows that the con- 
tribution from the term proportional to plasma pressure is 
the largest for both the upstream and downstream bound- 
aries. During the early stages of reconnection the Poynting 
flux out of the layer is comparable to the kinetic energy ef- 
flux, but as reconnection continues to develop the Poynting 
flux becomes small (<15%) compared to the kinetic energy 
efflux. Evaluating the outflow component of the conserva- 
tion of momentum integral given in equation (8) again shows 
that the plasma pressure term is dominant. Early in time, 
the downstream magnetic pressure due to the vertical mag- 
netic field is comparable to the momentum flux out of the 
layer but becomes small in comparison as reconnection de- 
velops and the outflow velocities increase with time. Mag- 
netic tension forces associated with the upstream boundaries 
are of the same order as the momentum flux exiting each 
side of the layer but are smaller than the contribution from 
terms proportional to pressure. The tension forces towards 
each downstream region are symmetric to within ~5-30% 
for case A, but for case B, the tension force directed to- 
wards the obstructing wall is ~2-3 times larger than the 
tension force directed towards positive x. The full evalua- 
tion of these surface integrals shows that equations (18) and 
(20) representing conservation of mass and energy can be 
used to successfully describe the scaling of steady magnetic 
reconnection with asymmetry in the outflow direction. For 
modest aspect ratio current sheets such as those associated 
with the Earth's magnetotail or spheromak merging, contri- 
butions to tension along the boundary can be important for 
momentum balance in the outflow direction and should be 
considered further in future work. 

5. Summary and Conclusions 

Magnetic reconnection with asymmetry in the outflow di- 
rection occurs in many systems in nature and in the labo- 
ratory, including planetary magnetotails, coronal mass ejec- 
tions, flux cancellation events, laboratory reconnection ex- 
periments, astrophysical disks, and magnetized turbulence. 
In this paper, we perform a control volume analysis to de- 
scribe long and thin current sheets with asymmetric down- 
stream pressure and test these scalings using resistive MHD 
simulations of driven reconnection. 

In section 3, we derive a set of scaling relationships which 
describe steady-state magnetic reconnection in a current 
sheet with asymmetry in the outflow direction without ex- 
plicitly specifying the dissipation mechanism. We derive ex- 
pressions for the outflow velocity for both the compressible 
and incompressible cases that do not directly depend on 
the dissipation mechanism. When resistive dissipation is as- 
sumed, we present an expression for the reconnection rate 
that depends on the outflow velocities from both sides of 
the current sheet. Together, these relations show how the 
outflow velocities and reconnection rate depend on a com- 
bination of upstream and downstream parameters. In the 
presence of asymmetric downstream pressure, it is possi- 
ble to have one Alfvenic jet and one sub- Alfvenic jet rather 
than two bidirectional Alfvenic jets. The reconnection rate 
is greatly reduced only when outflow from both sides of the 
current sheet is blocked. This helps explain results by Oka 
et al. [2008], who find that the presence of an obstacle on 
one downstream side of the current sheet does not greatly 
impact the reconnection rate. 

In a steady state, the magnetic field null and flow stagna- 
tion point overlap only in the absence of pressure gradient 
forces at the magnetic field null. When there is a pressure 
gradient, the magnetic field null is located on the side of 



the flow stagnation point which allows magnetic tension to 
counter the non-electromagnetic forces at the flow stagna- 
tion point. The position of the flow stagnation point can 
be estimated using conservation of mass when the upstream 
density and inflow velocity are approximately uniform. The 
position of the magnetic field null relative to the flow stagna- 
tion point can be estimated using a Taylor expansion around 
the flow stagnation point. 

To test the scaling relations derived in this paper, we per- 
form two-dimensional resistive MHD simulations of driven 
reconnection using the setup of MRX in linear geometry. 
Asymmetry in the outflow direction develops because one 
downstream wall is closer to the current sheet than the 
other downstream wall. The driving mechanism of MRX 
constrains the current sheet position between the flux cores 
and limits current sheet motion. Data extracted from this 
test show good correspondence with the scaling relations ap- 
proximating conservation of mass, momentum, and energy. 
The solution of equation (21) for outflow velocities shows 
reasonable agreement but increased scatter since the roots 
of high order polynomials can be sensitive to small errors 
in the coefficients. The reconnection electric field strength 
and the flow stagnation point position are well predicted 
by equations (32) and (37). The position of the magnetic 
field null is not well predicted by equation (45) due to the 
presence of a local pressure maximum near these two points. 
Exact evaluation of the integrals show that most of the as- 
sumptions of the model are met but that there is a non- 
negligible contribution due to tension along the boundary 
of the current sheet. 

Laboratory plasma experiments such as MRX, SSX, and 
TS-3/4 provide an excellent opportunity to study the impact 
of asymmetry on the reconnection process. The pull mode of 
operation in MRX is well-suited to investigate reconnection 
with asymmetry in the inflow direction due to cylindrical ge- 
ometry effects [Murphy and Sovinec, 2008]. However, effects 
related to downstream pressure may need to be incorporated 
into the scaling relations of Cassak and Shay [2007]. SSX, 
TS-3/4, and the push mode of operation in MRX can be 
used to study the impact of asymmetry in the outflow di- 
rection. The effects of asymmetry in the outflow direction 
(including current sheet motion) can be further studied by in 
situ measurements in the magnetotail and by observations 
of solar reconnection phenomena such as flux cancellation 
events, chromospheric jets, solar flares, and coronal mass 
ejections. 

The model developed in this paper assumes steady-state 
two-dimensional antiparallel reconnection in a high aspect 
ratio current sheet. Refinements or alternatives to this anal- 
ysis would benefit from the inclusion of time-dependent ef- 
fects such as current sheet motion and plasmoid formation. 
Of particular interest are what determines the rate of X-line 
retreat as seen in simulations by Oka et al. [2008] and how 
the current sheet structure and dynamics are changed due 
to current sheet motion. Three-dimensional effects have the 
potential to enhance the ability of plasma to exit the cur- 
rent sheet [e.g., Lazarian and Vishniac, 1999; Sullivan and 
Rogers, 2008; Shimizu et al, 2009]. Future analyses should 
consider the uneven contribution of magnetic tension for 
modest aspect ratio current sheets. This would allow the ef- 
fect noted by Galsgaard et al. [2000] , Galsgaard and Roussev 
[2002] , and Murphy and Sovinec [2008] , in which asymmetric 
outflow develops because the X-point is displaced towards 
one end of the current sheet, to be quantified. 
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